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Abstract: The model of three-body Borromean halo nuclei breakup was described by using standard phase 
space distributions and the Monte Carlo simulation method was established to resolve the detection problem 
of two neutrons produced from breakup reaction on the neutron wall detector. For ®He case, overall resolution 
(JBj, for the Gaussian part of the detector response and the detection efficiency including solid angle acceptance 
with regard to the excitation energy Ek are obtained by the simulation of two neutrons from ®He breakup into 
the neutron wall. The effects of the algorithm on the angular and energy correlations of the fragments are 
briefly discussed. 
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1 Introduction 

The neutron wall [l|, 0], which is located at the 
external target experimental hall of CSRm at the In- 
stitute of Modern Physics, Chinese Academy of Sci- 
ences, has been built to measure neutrons in near- 
relativistic heavy-ion collisions by the time of flight 
(TOF) method. The main aim of the neutron wall 
is the complete measurement of neutron momenta to 
reconstruct the invariant mass of the nucleus, which 
is feasible due to strong kinematical focussing in the 
case of neutron emission from excited projectiles. 

There is a specific problem which is related to the 
neutron detection. A number of the neutron wall's 
submodules are fired by neutrons impinging on it and 
a pattern recognition algorithm must be employed in 
order to disentangle multiple neutron hits. The al- 
gorithm and its performance of one neutron, under 
the circumstances of an simulation very similar to 
the present one, is described in Ref. 3]. When two 
neutrons need to be detected by the neutron wall at 
the same time, the effects of the algorithm on the an- 
gular and energy correlations of the fragments must 



be studied in detail, and the solid angle acceptance of 
two neutrons also need to be considered. The main 
effect takes place in a reduced double-hit recognition 
capability in the case where two neutrons interact in 
close vicinity to each other in the neutron wall. In 
three-body Borromean halo nuclei Q breakup exper- 
iments, it has to be corrected for such detection de- 
ficiencies on the basis of realistic event simulations, 
adjusted to the experiments 

One simple three-body Borromean halo nuclei 
breakup model was presented in the following. As 
an example, ®He breakup reaction will be studied to 
confirm the simulation method and the recognition 
algorithm in this paper. 



The model of three-body 
romean halo nuclei breakup 



Bor- 



Each event of three-body Borromean halo nuclei 
breakup is characterized by nine parameters: three 
components of the momentum vectors (ki,k2,k3) for 
the two neutrons and the core-particle in the projec- 
tile rest frame 0j. In a three-body continuum state, 
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the fragment relative motion can be described by two 
relative Jacobi momenta in the following: fik^ be- 
tween two constituents and Kky between the third 
constituent and the center of mass (cm.) of the first 
pair. Then, the normalized Jacobi momentum coor- 
dinates are constructed as 



— i-ixi )j 
mi 1712 
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(mi -f ma) ' 



ki-f ka 



ma 



_ m,3(mi-|-m2) , , 



kc.m. — ki + ka + ka. 

where and ?iki, i = 1, 2, 3 are masses and momen- 
tum vectors, respectively, of the particle i in the pro- 
jectile rest frame and fi^, fJ-y are the reduced masses, 
M = 7ni+ m2 -|- rris . 

The sum of kinetic energies for the fragments' rel- 
ative motions, = E.^ + Ey = ?i^k|/2^^ + ff]iy/2fiy 
is equal to the excitation energy Ek of the three- 
body Borromean halo system above threshold. The 
energy Ek fixes only a total phase volume which is 
accessible for the fragments and fragment kinetic en- 
ergies have continuous distributions within this vol- 
ume. Two variables need to be introduced for future 
reference: the variable e — Ex/E^ [Ey/E^ = 1— e) 
describes the share of the relative kinetic energy re- 
siding within a pair of particles at energy E^ , and the 
variable cos^^^^ = (k^-kj^) describes fragment angular 
correlations. In the cm. coordinate system, where 
the total nucleus is at rest (kc.m. = 0), the momen- 
tum of the third particle h\<i-j, is equal to ftky. 



O"- ! O 




T-system 



Fig. 1. The T and Y Jacobi coordinate sys- 
tems for three-body Borromean halo nuclei 
with two identical nucleons and the core C. 

As is shown in Fig. [H the case where two of the 
three constituents are identical particles is discussed , 
implying only two different Jacobi coordinate system. 
The T coordinate system is 'cluster' representation. 



where the two neutrons have index 1 and 2 (coor- 
dinate kj,) and the core 3 is connected to the cm. 
of the two neutrons, 1 and 2 (coordinate k^,). The 
Y coordinate system is 'shell-model' representation, 
where the core have index 1 and is coupled to the 
neutron with index 2 (coordinate k^). The second 
neutron 3 is then connected with the cm. of 1 and 2 
(coordinate kj^). 

For the event simulation, it is assumed that the 
available kinetic energy is distributed among the core 
particle and the two neutrons according to standard 
phase space distributions and zero momentum trans- 
fer to the excited projectile. The fractional energy 
and angular correlations in the core + n+n system at 
fixed Ek can be described by a probability distribu- 
tion W(e,6'). It represents the probability density of 
finding the whole system in a final state with a con- 
figuration having fixed values of e and 6 in the coor- 
dinate system T and Y. The probability distribution 
is normalized to unity. The distributions about the 
kinetic energy shared in the subsystems and the polar 
angle between Jacobi vector k^, and kj, in the model 
are given by 



dW/de 



dW/dcosOx 



v/e(l-e), 



constant. 



(2) 



3 The Monte Carlo simulation 
method 

In order to obtain the experimental performance 
of the neutron walL a Monte Carlo simulation code 
based on Geant4 [7| was implemented. "^He is one of 
simple three-body Borromean halo nuclei with a lot of 
experimental results 0, Q which can be a good refer- 
ence for the simulation code. Therefore, the response 
of the neutron wall with selected excitation energies 
of ''He, disintegrating into a-|-n-f n, was studied by 
simulation. It can be prepared for the other similar 
three-body Borromean halo nuclei experiments on the 
neutron wall in the future. 

The relative setup in the simulation is the same 
to the one described in Ref. Q but the distribu- 
tion of the neutrons that incident into the neutron 
wall. The neutrons are generated from the model 
mentioned above without considering the simulation 
of the a particle. It can be considered that the a 
particles are completely detected with no deviation. 
The three-body phase space Ex Ey dEx dEy and E^ 
are invariants, i.e. independent of the Jacobi system. 
Therefore, T Jacobi coordinate system for ''He with 
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the two neutrons and the a particle is selected for 
the simulation. The selected excitation energies of 
''He breakup into a + n + n are 0.5, 1, 1.5, 2, 3, 4, 5 
and 6 MeV according to Ref . [8] . The energy of '^He 
ion is 300 MeV/nucleon. 

4 Results and discussion 

Taking into account intrinsic detection efficiencies, 
position resolutions, the TOF resolution, and finite 
acceptances for the two neutrons, the response for a 
given excitation energy of ^He and subsequent decay 
into a+n+n is derived by event simulation described 
above. The algorithm to disentangle the two neu- 
trons impinging on the neutron wall is identical for 
the analysis of simulated events in Ref. 




Fig. 2. (color online) (Upper panel): Resolu- 
tion CT_Bj, for the Gaussian part of the neutron 
wall response (right hand scale, square sym- 
bols) and detection efficiency e{Ek) including 
solid angle acceptance (left hand scale, cir- 
cle symbols). (Lower panel): Response of the 
neutron wall for selected excitation energies 
Ek = 0.5, 1, 1.5, 2, 3, 4 and 5 MeV of ^He, 
disintegrating into a + n + n (see text). 

With this simulation procedure, we obtain the de- 
tection efficiencies shown in the upper panel of Fig. [2] 
It is observed that two neutrons impinging on the de- 
tector are well recognized in about 30% of the events 
on the maximum. We also like to point out, that even 
two neutrons incident with zero relative angle can be 



disentangled in about 15% of such cases. This can 
be qualitatively understood by comparing the depth 
of the " shower" volume (about 60 cm) with the total 
depth of the detector (100 cm) and taking into ac- 
count the mean interaction length of about 54 cm for 
energetic neutrons in the neutron wall. Therefore, 
a differential efficiency dependent on the neutron- 
neutron relative angle arises which e.g. is reflected, 
and to most extent determines the efficiency shown 
in Fig. [2l The apparent decrease in efficiency, at low 
is due to the limited capability to resolve two 
neutrons with a small relative distance in the neu- 
tron wall. The decrease at higher excitation energies 
is due to the finite solid angle acceptance. 

The same procedure also delivers response func- 
tions which can be well described by Gaussian distri- 
butions (Fig. [21 lower panel). Small non-Gaussian 
wings are presented on a few percent level which 
can be neglected under most circumstances. For the 
Gaussian part, we observe a variance aE^iMeV) « 

0. 13(£;fe)^ °^ with Ek in MeV as shown in the upper 
panel of Fig. [5] 

A comparison of the angular and energy corre- 
lation distributions for the excitation energies Et^ = 

1, 2 and 5 MeV obtained in our analysis with those 
from the model described above is shown in Fig. [3l 
The left column shows the angular correlation distri- 
butions depending on cosO^y in the T system; the 
right column gives the energy correlation ones. For 
the angular correlations, the difference between the 
detected and reconstructed event distributions after 
the analysis, and the model ones at the three energies 
is not obvious, but the most striking feature is that 
the reconstructed results in a concave angular distri- 
bution for Ek= 1 MeV (Fig. [Sfthe top panel in the 
left column)). This is the result of the low excitation 
energy, which is relative to the reconstructed error. 

For the energy correlation distributions, shown in 
the right column, the decrease in the probability dis- 
tributions between the model and the detected events 
at low £ for = 1 and 2 MeV, is due to the limited 
capability to resolve two neutrons with a small rela- 
tive distance in the neutron wall. And, the decrease 
at high e for Ek= 5 MeV, is due to the finite solid an- 
gle acceptance. There is almost no difference between 
the detected and reconstructed event distributions at 
the three energies. It demonstrates that the algo- 
rithm to disentangle the two neutrons impinging on 
the neutron wall is correct on some level. 
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Fig. 3. (color online) Projections of the probability distribution W{e,6) in T Jacobi coordinate system ob- 
tained in the selected excitation energies Ek = 1, 2 and 5 MeV. The black solid lines correspond to the 
probability distributions of the standard phase model, see text; the red dashed lines correspond to the ones 
of the detected events; the full triangles correspond to the ones of the reconstructed events after the analysis. 



5 Conclusion 

The model and Monte Carlo simulation method 
for identifying multiple neutron hits were established 
to resolve the detection problem of two neutrons 
in three-body Borromean halo nuclei breakup reac- 
tion. The three-body breakup ^'He — !■ a + n + n 
was studied as an example by simulation, using 300 
MeV/nuclueon ®He ion according to the model's dis- 
tributions with selected excitation energies. In the 
simulation, the response of the neutron wall was ob- 
tained by the two neutrons based on the three-body 



Borromean halo nuclei breakup model. We discussed 
the distributions of the efficiency t{Ek) and resolution 
CTbj. with regard to the excitation energy of ''He, 
obtained from the event simulations. The agreement 
between the detected event distributions and the re- 
constructed ones after the analysis is reasonably good 
at the three energies. Therefore, the algorithm to dis- 
entangle two neutrons impinging on the neutron wall 
is reliable, which could be helpful for the experiments 
with the detection of neutrons at the external target 
experimental hall in the future. 
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